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Abstract. We present the first 3D Monte Carlo photoionisation code to include a 
fully self-consistent treatment of dust radiative transfer (RT) within the photoionised 
region. This is the latest development of the recently published pure photoionisation 
code MOCASSIN (Ercolano et al., 2003a) and it is currently undergoing several bench- 
marking tests. The preliminary results of these tests are presented in these conference 
proceedings. The new code provides the ideal tool for the analysis of dusty ionised 
regions showing asymmetries and/or density and chemical inhomogeneities 



1 Introduction 

Recent advances in instrumentation allow the geometry of many PNe (including 
those in the SMC and LMC) to be resolved. It is now clear that the assump- 
tion of spherical symmetry is not justified in most cases. This imposes serious 
limitations on the application of ID codes, with the risk of obtaining incorrect 
abundance determinations, resulting in misleading conclusions. The 3D pho- 
toionisation code MOCASSIN (MOnte CArlo Simulations of Ionised Nebulae, Er- 
colano et al. 2003a) was designed to remedy these serious shortcomings and has 
been successfully applied to the study of several PNe (Ercolano et al. 2003b, c; 
Ercolano et al. 2004). 

The presence of dust grains in ionised environments can have significant 
effects on the radiative transfer, as the grains compete with the gas for the 
absorption of UV photons, as well as being heated by nebular resonance line 
photons. These processes can only be treated properly by incorporating the 
scattering, absorption and emission of radiation by dust particles that are mixed 
with the gas in the photoionised region. Moreover, the accurate determination of 
dust temperatures and spectral energy distributions (SEDs) can only be achieved 
by treating discrete grain sizes and different species separately. All of the above 
has already been implemented in the new version of MOCASSIN, which is currently 
undergoing benchmarking tests, as described below. 



2 Benchmarking Tests 



The dust and gas version of MOCASSIN is benchmarked against the popular ID 
code, dusty, for the spherically symmetric, homogeneous test models described 
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by Ivezic et al. (1997). The input parameters for these pure-silicate models are 
given in their paper and the grain optical constants used are those included in 
DUSTY's library, calculated from the data of Draine & Lee (1984). 

Figure 1 shows the comparison between MOCASSIN - crosses - and DUSTY 
(Ivezic et al., 1997) - dashed line- for the p=0 (see Ivezic et al., 1997) and 
Tljum = 10. The emerging SED in relative flux units are shown on the left panel 
and the radial grain temperature distribution on the right panel. Taking into 
account that a completely different approach is used by the two codes, the results 
obtained by the benchmarking are in good agreement. 2D disk benchmarks are 
currently being carried out for the pure dust case, whilst a dust and gas model 
of the PN NGC 3918 is also being constructed (Ercolano et al. in preparation). 

While a number of Monte Carlo dust radiative transfer codes already exist 
(e.g. Bjorkman & Wood, 2001), MOCASSIN is currently the only 3D code capable 
of self-consistently treating dust RT within the photoionised region, providing 
the ideal tool for the modelling of non-spherical, density and/or chemically in- 
homogeneous dusty photoionised environments. 
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Fig. 1. Spherically symmetric, homogeneous benchmark model. See text for details. 
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